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ABSTRACT 

Weak lensing by large-scale structure provides a unique method to directly measure matter fluctuations 
in the universe, and has recently been detected from the ground. Here, we report the first detection 
of this 'cosmic shear' based on space-based images. The detection was derived from the Hubble Space 
Telescope (HST) Survey Strip (or "Groth Strip"), a 4' x 42' set of 28 contiguous WFPC2 pointings with 
/ < 27. The small size of the HST Point-Spread Function (PSF) affords both a lower statistical noise, and 
a much weaker sensitivity to systematic effects, a crucial limiting factor of cosmic shear measurements. 
Our method and treatment of systematic effects were discussed in an earlier paper (Rhodes, Refregier 
& Groth 2000). We measure an rms shear of 1.8% on the WFPC2 chip scale (1.27'), in agreement with 
the predictions of cluster-normalized CDM models. Using a Maximum Likelihood (ML) analysis, we 
show that our detection is significant at the 99.5% confidence level (CL), and measure the normalization 
of the matter power spectrum to be agf^ 48 = 0.51^q 17 > m a ACDM universe. These 68% CL errors 
include (Gaussian) cosmic variance, systematic effects and the uncertainty in the redshift distribution 
of the background galaxies. Our result is consistent with earlier lensing measurements from the ground, 
and with the normalization derived from cluster abundance. We discuss how our measurement can be 
improved with the analysis of a large number of independent WFPC2 fields. 

Subject headings: cosmology: observations - gravitational lensing - cosmology: large-scale structure of 
universe - cosmology: dark matter 
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1. INTRODUCTION 

Weak gravitational lensing has now been established as 
a unique technique to measure the mass distribution in the 
universe (see Mellier 1999; Bartelmann & Schneider 2000 
for reviews). Recently, several groups detected this effect 
in the field by measuring the coherent distortions it in- 
duces in the shape of background galaxies (Wittman et al. 
2000; van Waerbeke et al. 2000; Bacon, Refregier & Ellis 
2000a (BRE); Kaiser, Wilson & Luppino 2000; Maoli et 
al. 2000) . This opens wide prospects for the measurement 
of the distribution of dark matter on cosmological scales 
and of cosmological parameters. 

In this paper, we present the first detection of this 'cos- 
mic shear' based on space-based images. For this purpose, 
we use the HST Survey Strip (or "Groth Strip" ), a 4' x 42' 
survey consisting of 28 contiguous WFCP2 fields with a 
limiting magnitude of I « 27 (Groth et al. 1994). A 
shear measurement method adapted to HST images and 
a detailed study of systematic effects for the Strip was 
presented in an earlier paper Rhodes, Refregier & Groth 
(2000, RRG; see also Rhodes 1999). The main advantage 
of HST images for weak lensing is the small PSF (0.1" 
compared to ~ 1" from the ground), which allows a larger 
surface density of resolved galaxies and makes the shear 
measurement much less sensitive to the PSF smearing. 
We apply the results to the Strip and search for a lens- 
ing signal using a Maximum Likelihood (ML) approach. 
We then derive constraints on the amplitude of the mass 
power spectrum and compare it to other measurements. 

The paper is organized as follows. In ^ we summarize 
the main theoretical results which serve as a reference for 



our measurement. In ^ we describe the relevant proper- 
ties of our data set. Our method for measuring the shear 
and our treatment of systematic effects are described in §[l| 
Our results are presented in §|| and discussed and summa- 
rized in SpL 



2. THEORY 

We measure the average shear 7 y i within each WFPC2 
chip, which are square cells of side a = 1.27'. The angu- 
lar fluctuations of the resulting smoothed shear field are 
characterized by the shear correlation functions Ci(0) = 
(7-' (0)t£(#)}, where i = 1,2 and the brackets refer to an 
average over pair of cells separated by an angle = \9\. 
The shear components 7J are measured in a rotated coor- 
dinate system whose x-axis is aligned with the separation 
vector between the two points (see e.g. Heavens, Re- 
fregier & Heymans 2000). These correlation functions can 
be computed for any cosmological model using (see eg. 
BRE) 
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[J (W)± MW)], 1-1,2, (1) 



where Wi is the window function of a square cell. The 
shear power spectrum Ci is defined as in BRE and can 
be evaluated from the evolution of the non-linear matter 
power spectrum using the fitting formula from Peacock & 
Dodds (1997). In this paper, we will consider a ACDM 
model with f2 m = 0.3, f^A = 0.7, V — 0.25, and consider 
several values of as, the normalization of the mass fluctu- 
ations on 8 ft -1 Mpc scales. 



2 



The variance of the cell-averaged shear is of cns = 
(|7| 2 ) = Ci(0) + C 2 (0). For the ACDM model, the rms 
shear in a cell is 

/(Xa x 1.24 /Q \ - 60 , z N 0.85 

^ 0.0197 (2) («=) (£) , (2, 

where z m is the median redshift of the background galax- 
ies, whose redshift distribution was taken to be of the form 
of Equation (0) below. The scaling relations in Equa- 
tion (0) provide an excellent approximation in our range 
of interest, namely, 0.6 <; cr s < 1.4 and 0.7 < z m < 1.2. 

3. DATA 

The Hubble Space Telescope Survey Strip (sometimes 
referred to as the 'Groth Strip', or the 'Groth-Westphal 
Strip') was taken during March and April 1994 using the 
Guaranteed Time Observations of one of us (EJG) and 
the WFPC-1 Instrumentation Definition Team (Groth et 
al. 1994). The Strip was observed in two passbands, V 
(F606W) and I (F814W). For 27 of the pointings in the 
Strip, total exposure times were 2800 and 4400 seconds in 
the V and I bands, respectively. One of the Strip point- 
ings (the Strip Deep Field) had exposure times of 24,400 
and 25,200 seconds, arranged in 4 dithers. We used only 
the I-band images to study weak lensing. The Strip im- 
ages consist of 4 individual exposures with equal exposure 
times. Here, we chose to use only one of the four dithers 
for the Deep Field, thus avoiding the complications in the 
PSF correction that dithering would introduce (see RRG). 

We used the Faint Object Classification and Analysis 
System (FOCAS; Jarvis and Tyson, 1981) within IRAF 
to create a catalog of objects in the Strip. We determine 
the sky background and its standard deviation using the 
IRAF task Imarith.. This was found to be more reliable 
than the corresponding algorithm within FOCAS. Mag- 
nitude zero points are those defined by Holtzman et al. 
(1995). A FOCAS catalog was created for each filter, and 
a matched catalog was created for objects that had I and 
V positions that coincided within a radius of about 5 pix- 
els. There were approximately 10,000 matched objects 
with I< 26. Of these objects, approximately 4000 were 
classified as galaxies, had well defined moments and were 
sufficiently large to be used in our weak lensing study. 
We were careful not to double count galaxies that fall in 
the small overlap region between consecutive fields in the 
Strip. 

Koo et al. (1996) have shown that a small subsample 
(24 galaxies) of Strip objects matches an extrapolation of 
the CFRS sample redshift distribution. More recent red- 
shift measurements of galaxies in the Strip by the DEEP 
collaboration (DEEP Collaboration, 1999) agree with the 
preliminary findings in Koo et al. (1996). Assuming the 
functional form of the redshift distribution found in the 
CFRS, 

^ oc *V<«/*>" (3) 
dz 

and using the median redshift of z m ~ z a — 0.5 for ob- 
jects in the range 17 < I < 22 found in that survey (Lilly 
et al. 1995), we calculate a median redshift for objects 
in the Strip. The galaxy sample we use to study weak 



lensing (/ < 26 and size d > 1.5 pixels) has a median 
magnitude of I — 23.6. Extrapolating from the CFRS 
sample using the Strip number counts, this corresponds to 
a median redshift of z m = 0.9 (Rhodes 1999, Groth and 
Rhodes, 2001). Photometric redshifts for galaxies in the 
Hubble Deep Fields (HDF-North and HDF-South) have 
been measured by Lanzetta et al. (1996). A linear fit to 
the median redshift versus median magnitude plot yields 
z m = 0.87 for the HDF-South and z m = 1.00 for the HDF- 
North for a median magnitude of / = 23.6, where we have 
used the magnitude conversion I = —0.49 + Iab- This is 
consistent with what we found above using an extrapola- 
tion of the CFRS sample. Therefore, we adopt a median 
redshift of z m = 0.9 ± 0.1 for the lensing analysis in this 
paper. 

4. PROCEDURE AND SYSTEMATIC EFFECTS 

The procedure we use for extracting galaxy ellipticities 
and shear from the source images is described in detail in 
RRG (see also Rhodes 1999). This method is based on 
one introduced by Kaiser, Squires, and Broadhurst (KSB, 
1995), but modified and tested for application to HST im- 
ages. The KSB method was shown to be adequate for 
current ground based surveys (Bacon, et al. 2000b; Erben 
et al. 2000). In RRG, we used numerical simulations and 
an array of tests to show that our method was sufficiently 
accurate for WFPC2 images. 

We use our method to correct for two effects: camera 
distortion and convolution with the anisotropic PSF. Cor- 
rections are done using moments measured with a Gaus- 
sian weight function whose size in pixels, oj, depends on 
the size of the object as ui — max(2, tJA/tt), where A is 
the area of the object in pixels. The minimum size of 
two pixels was found to be the optimal weight function 
size for stellar objects using both actual WFPC2 data and 
simulations using artificial stellar images created using the 
program Tiny Tim (Krist and Hook 1997). 

Camera distortions are corrected using a map derived 
from stellar astrometric shifts (Holtzman, et al., 1995). 
PSF corrections are determined from HST observations 
of two globular clusters (M4 and NGC6572). Finally, we 
measure the ellipticities (e,) of the galaxies in the Strip 
and convert them into shear estimates using 7$ = G~ l ti, 
where G is the shear susceptibility factor given by Eq. (30) 
in RRG. 

RRG have shown that the residual systematic effects can 
be greatly minimized if (1) the shear is averaged over entire 
chips, and (2) only "large" galaxies with a (convolved) rms 
radius greater than 1.5 pixels (or about 0.15") are selected. 
Since the PSF anisotropy is mainly tangential about the 
chip centers, condition (1) ensures that the mean PSF el- 
lipticity is small (e* ~ 0.02). Since the shapes of larger 
galaxies are less affected by the PSF, condition (2) en- 
sures that the impact of the PSF anisotropy on galaxy 
ellipticities e 9 is small. Specifically, RRG show that the 
ellipticity of the selected galaxies e g induced by the PSF 
ellipticity is reduced by a factor / rec j = e 9 /e* ~ 0.13. 

As shown in RRG, the residual systematics for this 
galaxy subsample are dominated by two effects. The first 
arises from the imperfect accuracy of the ellipticity cor- 
rections and shear measurement method. This effect re- 
sults in an uncertainty in the normalization of the shear 
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of about 4% and is easily incorporated at the end of our 
analysis. The second results from the time variability of 
the PSF. To quantify this effect, we measure the mean el- 
lipticity of the PSF in each chip for 4 stellar fields (M4 
observed at 2 different epochs, NGC6572, and a stellar 
field from the WFPC2 parallel archive). We find that the 
chip-averaged PSF ellipticity e* c varies by about 0.01 (rms) 
in chips c = 2 and 3 and 0.02 in c = 4. The variations 
appear to be stochastic and to be unrelated to the focus 
position of the telescope. These time variations must be 
corrected statistically since the small number of stars in 
the Strip precludes measuring the PSF in each field in- 
dividually. For this purpose, we measure the covariance 
C* c j C ' = cov[e* c , e* c ,] of the stellar ellipticities by averaging 
over all stellar fields, and convert this into the systemat- 
ics shear correlation matrix C S J, S - C , ~ f^G~ 2 C* c JC , which 
is subtracted from the galaxy shear correlation functions. 
The diagonal elements are the shear variance produced by 
the time variability of the PSF cr s 2 ys = (C*^ s lc + CZ% c )c - 
0.001 1 2 , where the average is over all chips. 



5. RESULTS 

To study the lensing statistics in the Strip, we first com- 
pute the mean shear j ic in each chip c, by averaging over 
all selected galaxies in the chip. This mean shear is the 
sum of contributions from lensing, noise and systematic 
effects, i.e. j lc = + 7?°™ + 7^ . From the dis- 

tributions of the galaxy shears within each chip, we can 
measure the covariance matrix Cf°J^, = cov[t"" c oisc , 7"°/ se ] 
of the mean shear components. Since the chips do not 
overlap, Cf® 1 ^, vanishes for c ^ c'. The diagonal ele- 
ments Cf°f° 



°noisc ic are the errors in the mean for 



the mean shear components j ic . The average noise vari- 

} c ~ 0.047 2 , where the 



ance is err 



' noise, 2c/c 



average is over all chips. The off-diagonal entries in the 
noise matrix are an order of magnitude below the diagonal 
ones and thus do not have much impact on the analysis. 
We nevertheless keep them for completeness. 

The shear pattern for the Strip is shown in Figure [lj 
Since the expected signal on the chip scale is about a\ ens — 
0.02 (see Eq.[2]), the signal-to-noise ratio ai cns /a no i se for 
each chip is about 0.4. The lensing signal must there- 
fore be searched for statistically. For this purpose, we 
first consider the mean shear (7fc) c averaged over all 
chips in the Strip. We find (7!) = 0.005 ± 0.004 and 
(72) = 0.005 ± 0.004, where the errors are la errors in 
the mean. This is consistent with zero, as expected given 
the length of the Strip. This nonetheless confirms that we 
are not subject to systematic effects on large scales. 

1— ilQV. Ellipticity 
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Fig. 1. — The ellipticity of each field and chip in the Survey 
Strip. The Strip measures about 4 by 42 arcminutes. 

We next consider the variance of the shear on the 
chip scale. In general, the total observed shear variance 
tr 2 ot = (|7j C | 2 ) is the sum of the contribution from lens- 
ing, from noise, and from systematic effects, i.e. of ot = 



CT iens + cr noisc+ 'sys- An estimator for the lensing variance is 
thus af c ns ~ cr t 2 ot - cr 2 oisc - cr s 2 ys (for more details, see BRE) . 
Using our estimate of a sys (see §||), we find a\ cns ~ 0.018, 
which is close to the expected value of 0.0197 in the ACDM 
model (see Eq. g]). Estimating errors on the shear vari- 
ance is complicated since the chips are contiguous and thus 
not independent. Instead of attempting to do this, we use 
a maximum likelihood approach to extract the lensing sig- 
nal on all scales simultaneously. 

To do so, we assume that the shear field obeys Gaussian 
statistics. Since we are averaging over a large number of 
galaxies (about 50) in each field, the central value theorem 
ensures that this assumption is valid for the noise contri- 
bution. It is difficult to establish whether the systematics 
are Gaussian; however, since they have a much smaller 
amplitude than the noise, it is sufficient to assume that 
they are so. The lensing shear field is known not to be 
Gaussian, especially on scales smaller than 10' which are 
dominated by the non-linear evolution of structures (Jain 
& Seljak 1997). However, as we have shown above, the 
noise dominates over the lensing signal in each individual 
chip. Therefore our Gaussian approximation should not 
underestimate our error bars by a large amount. 

With this approximation, the probability of measuring 
the mean shears j ic given a cosmological model M is 



P(l ic \M) 



y v c- 



{2n) N \C\i 



(4) 



where N = 84 is the number of chips, and g — (7 11: 7 2 i, ■ ■ 
•,7 1Ar ,7 2 jv) is the data vector composed of the 2N mea- 
sured shears. The total correlation matrix C is given by 
C = C lens + C noiso + C sys . The systematics correlation 
matrix C sys is computed as explained in §[|, assuming 
that each field is independent as indicated by the stochas- 
tic nature of the PSF variations. The lensing correlation 
matrix (7 lcns depends on the cosmological model and can 
be computed by unrotating the rotated correlation func- 
tions C\(9) and 6*2(0) (Eq. m). Following Bayes' postulate, 
the likelihood L of the model M given the data is then 
L cx P(y ic \M), so that its logarithm is given by 



21nL ~ g T C~ 1 g + In \C\ + constant. 



(5) 



First, we test the null hypothesis, namely the absence 
of lensing corresponding to C lcns = 0. In this case, 
C = C noiso + C sys is a constant, and the log likelihood 
reduces to — 21nL = \ 2 + constant, where \ 2 = 9 T C~ x g. 
For the Strip, we find a reduced x 2 of 1.31 with 2 x 84 de- 
grees of freedom. The null hypothesis is thus ruled out at 
the 99.5% confidence level, thus establishing our detection 
of the lensing signal. Note that this null test does does 
not rely on the assumption that the lensing shear field is 
Gaussian, because we have set it to zero in this case. 

Now that we have established the presence of a lensing 
signal, we can compute the constraints our data impose on 
cosmological parameters. To do so, we compute the Maxi- 
mum Likelihood (ML) estimator for CT 8 (fi m /0.3) - 48 . This 
was done by maximizing InL (Eq. |H), which was com- 
puted with Ci cns matrices corresponding to different val- 
ues of this parameter. To compute the error in <j% we per- 
formed 1000 random realizations of the Strip with shears 
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drawn from a multivariate Gaussian distribution with a co- 
variance matrix given by C = C* lcns + C noisc + C sys (with 
C lons at the MLE value of a 8 (fl m /0-3) OA8 of the Strip). 
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Fig. 2. — Distribution of the ML normalization of the power 
spectrum asQ. 1 ^ 8 derived from 1000 random realizations of the 
Strip. The normalization corresponding to the real Strip is 
shown as the vertical dashed line. The 68% confidence level 
is shown as the shaded region. The full 68% confidence inter- 
val which includes (Gaussian) cosmic variance, the uncertainty 
in the median redshifts of the background galaxies and that 
arising from systematic effects is shown as a horizontal bar. 
It can be compared to the corresponding limits derived from 
the William Herschel Telescope (WHT) cosmic shear survey 
of BRE, and to that from cluster abundance (Pierpaoli et al. 
2000), shown as marked horizontal bars. 

The resulting distribution of the ML values is shown 



0.48 _ 



0.9(T 



in Figure |2j. We thus obtain cr 8 (f2 m /0.3) 
where the error bars correspond to the 68% confidence 
interval. After propagating errors induced by the uncer- 
tain median rcdshift z m (see §[|) of the galaxies and by 
the shear measurement method, our final 68% confidence 
range becomes <jg(tt m /0.3) OA8 = 0.901q3o anc ^ ^ s snown 
as a horizontal bar in Figure 0. This value is consistent 
with the William Herschel Telescope (WHT) cosmic shear 
survey of BRE, who found cr 8 (n m /0. 3)° A8 = 1.5 ± 0.5, 
as shown by the associated horizontal bar on the Figure. 
Our results are also consistent with the amplitude of the 



cosmic shear signals found by other groups (Wittman et 
al. 2000; van Waerbeke et al. 2000; Kaiser et al. 2000; 
Maoli et al. 2000). However, since these groups do not 
include cosmic variance and redshift uncertainty in their 
error bar estimates, their results can not be directly com- 
pared with ours. The normalization from cluster abun- 
dance cr 8 (rW0.3) - 60 = 1-O19±o:o76 (Pierpaoli et al. 2000) 
is consistent with ours and is also shown on the figure (af- 
ter ignoring the small difference in the exponents for tt m ). 



6. CONCLUSIONS 

We have detected, for the first time, cosmic shear us- 
ing space-based images. Our detection in the HST Survey 
Strip is significant at the 99.5% confidence level. Using 
a ML method to search for a lensing signal on all scales 
simultaneously, we derived a normalization of the matter 
power spectrum of crgO^ 48 = 0.51±°;^. This 68% con- 
fidence interval includes statistical noise, (Gaussian) cos- 
mic variance, uncertainty in the redshift distribution of 
the galaxies and in our shear measurement method. This 
is consistent with ground-based cosmic shear surveys and 
with the normalization derived from cluster abundance. 

Our result demonstrates the power of space-based im- 
ages to measure cosmic shear. The small PSF in HST 
images affords both a larger surface density of background 
galaxies and a much reduced impact of the systematics in- 
duced by the PSF anisotropy. The latter advantage will 
be even more important for future, wider surveys whose 
sensitivity will be comparable to the amplitude of the sys- 
tematics. 

Our current results are chiefly limited by the relatively 
small area of the Strip. The signal-to-noise ratio for the de- 
tection of cosmic shear in Nf independent WFPC2 fields, 
in which we cannot search for lensing on all scales simul- 
taneously as we have done here, is approximately S/N 

K„ s /CT no i S o) 2 (3A7)^ ^6.8(^/500)5 (see BRE). Our 
detection can thus be significantly improved in the future 
by an analysis of existing disjoint WFPC2 fields. 
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